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Objectives 

The principal objective of this unit is to introduce differential equations 
of the first-order and to discuss their soluuons and some methods of 
solving them, 


After working through this unit you should be able to: 
(i) explain the meanings of the following terms: 

solution set of a differential equation, 
family of solution curves, 
separation of variables method, 
integrating factor method, 
Euler's method of solving first order differential equations ; 

(ii) sketch a family of solution curves for a given simple first order 

differential equation ; 

(iii) solve a given first order differential equation by either the separation 
of variables method or the integrating factor method (in appropriate 
cases); 

(iv) having completed (iii), find a particular solution with a given initial 
condition ; 

(v) apply Euler’s method, given a first order differential equation, an 
initial condition and a specified step length. 


Note 


Before working through this correspondence text, make sure you have 
read the general introduction to the mathematics course in the Study 
Guide, as this explains the philosophy underlying the whole course. 
You should also be familiar with the section which explains how a text 
is constructed and the meanings attached to the stars and other symbols 
in the margin, as this will help you to find your way through the text. 
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Glossary 


Terms which are defined in the glossary are printed in CAPITALS. 


DEPENDENT VARIABLE See INDEPENDENT VARIABLE. 


DIFFERENTIAL 
EQUATION 


EULER'S METHOD 


FAMILY OF 
SOLUTION CURVES 


FORMULA METHOD 


FUNCTION OF ONE 
REAL VARIABLE 


FUNCTION OF TWO 
REAL VARIABLES 


INDEPENDENT 
VARIABLE 


INITIAL CONDITION 
FOR A FIRST 

ORDER DIFFERENTIAL 
EQUATION 


INTEGRATING 
FACTOR METHOD 


ORDER OF A 
DIFFERENTIAL 
EQUATION 


SEPARATION OF 


VARIABLES 
METHOD 


SOLUTION CURVE 


SOLUTION OF A 
DIFFERENTIAL 
EQUATION 


vi 


A DIFFERENTIAL EQUATION is an equation which 
relates an unknown real function (which we wish to 
find), its derived functions, and possibly other 
known functions. Such an equation can alterna- 
tively be expressed in terms of the images (variables) 
under the functions. (See the detailed discussion in 
the text.) 


EULER'S METHOD is a numerical method of solution of 
@ DIFFERENTIAL EQUATION. (See the detailed discus- 
sion in the text.) 


The FAMILY OF SOLUTION CURVES Of a DIFFERENTIAL 
EQUATION is a graphical representation of the 
SOLUTION SET, 


A FORMULA METHOD is a method of solution of a 
DIFFERENTIAL EQUATION which produces a formula 
for the solution. 


A FUNCTION OF ONE REAL VARIABLE is a function 
with domain a subset of R, and codomain R. 


A FUNCTION OF TWO REAL VARIABLES is a function 
with domain a subset of R x R, and codomain R. 


If Q:1-—+q, then 1, the domain variable of Q, is 
called the INDEPENDENT VARIABLE, and q is called the 
DEPENDENT VARIABLE. 


An INITIAL CONDITION FOR A FIRST ORDER DIF- 
FERENTIAL EQUATION is a condition which picks out 
a particular element of the SOLUTION sET. Geo- 
metrically, it represents a given point through which 
the graph of that particular element passes. 


The INTEGRATING FACTOR METHOD is a special 
method of solution of a first order DIFFERENTIAL 
EQUATION. (See the detailed discussion in the text.) 


The ORDER OF A DIFFERENTIAL EQUATION is the order 
of the highest derived function involved. 


The METHOD OF SEPARATION OF VARIABLES is a 
particular method of solution of a first order 
DIFFERENTIAL EQUATION. (See the detailed discussion 
in the text.) 


A SOLUTION CURVE Of a DIFFERENTIAL EQUATION is 
the graph of an element of the SOLUTION SET. 


A SOLUTION OF A DIFFERENTIAL EQUATION is an 
element of the SOLUTION SET. 
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SOLUTION SET OF 
‘A DIFFERENTIAL 
EQUATION 


STEP LENGTH 


The SOLUTION SET OF A DIFFERENTIAL EQUATION is the 
set of all real functions which satisfy the DiF- 
FERENTIAL EQUATION. 


The STEP LENGTH is the width of interval used to 
determine the next point on an estimated SOLUTION 
CURVE of a first order DIFFERENTIAL EQUATION. 


30 
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Most textbooks go well beyond the scope of this text. We have therefore 
chosen just two to give you some idea of the number of classical methods 
available, and to give a different viewpoint of some of the material 
presented in this text. They are 


H. Betz, P. B. Burcham and G. M. Ewing, Differential Equations with 
Applications (Harper and Row, 1964). 
Parts of Chapters 1-4 are relevant. 


W. Kaplan, Elements of Ordinary Differential Equations (Addison- 
Wesley, 1964). 
Parts of Chapters 1-3 are relevant. 


For the numerical approach, we suggest: 


N. Macon, Numerical Analysis (John Wiley, 1963). 

Chapter 10 gives an in ion of the more formal approach to the 
numerical methods and their accuracy. You will probably find this last 
reference rather hard at this stage. 


General Note 


In this unit we shall use, as appropriate and convenient, any of the calculus 
notations mentioned previously in the course. For instance, if Q is a real 
function, such that 


Q:t-—- Qt) (te R) 


and we write q = Q(t), then we write the derived function as Q’, when the 
“function” notation suits us, or we express the derivative as 


ry) = 49 
WW=T 


using the “function image” notation or the Leibniz notation. You should 
by now be reasonably familiar with the function notation as we have used 
it throughout the course, and also with the distinction between the func- 
tion Q and the variables t and q, where t is a general element of the domain 
of Q, and q is a general element of the image set. Our objective in using 
these various notations is to enable you not only to read this course, but 
also to read any other texts that you may come across now and in further 
studies. 
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24.1 INTRODUCTION 


24.1.0 What is a Differential Equation? 
In Unit /2, Differentiation I we saw that the derivative of a real function Q 
represents the rate of change of the image value under Q, or the slope of 
the (tangent to the) pictorial graph of Q; that is, Q’(1), (1 € R), is the rate of 
change of Q(r) at t. Given the function Q, we can derive Q' by using one of 
the standard rules. For example, if 

Q:-— et - 2° (te R), 
then 

Ost—>e' — 67 (te RR). 
This is straightforward. 


The problem is appreciably harder when we are given the derived function 
Q' and are asked to determine Q. This is the problem of integration and we 
have developed some tools to cope with this (Unit 9 and Unit /3, Integra- 
tion | and 11). Essentially, there are two main points. In the first place, the 
differentiation operator D is a function, but the reverse of D is nota 
function. We coped with this by introducing a constant function and 
obtained the images of Q’ under the reverse mapping in the form 


any primitive + constant function. 
For example, if we are given 

Ot—2t (teR) 
then 

Q:u-—SP +e (rER), 


where ¢ is a real number, and for each choice of c we get an image of Q’ 
under the integration mapping. The second important point is that very 
often we could not find a simple expression for the integral (or, even if we 
could, the labour involved was prohibitive). 


For example, if 


(teR*) 


Q:t— 


then 

Cn—+? 
In this unit we are going to make our integration process still more 
difficult ; not because it is “good for the soul”, but because of the wide field 
of application of the development. We are on the threshold of the subject 
known as differential equations, which has a considerable literature 
associated with it and on which much research is still undertaken today. 


A differential equation of the type which we consider in this unit is an 
equation involving an unknown real function and its derived functions, 
and possibly other known functions. Thus if f is a real function, 


L+Y=0 

(where O is the function x-— 0 (x € R)) is a differential equation, as is 
S’ + (xx?) x f = xin x. 

The simplest form of differential equation is an equation of the form 
ae 

where g is a known real function; its solution presents us with precisely the 

integration problem we have discussed previously. We shall discuss more 
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fully in section 24.1.2 what we mean by the solution of a differential equa- 
tion, but we note here that the solution of such an equation is a function (or 
set of functions). This contrasts with most previous equations we have 
solved in this course, where in each case the solution has usually been a 
number (or a set of numbers). In general, in mathematics an equation is a 
relation defined on a set of elements and can be considered as defining a 
subset consisting of those elements which belong to the solution set of the 
relation. (For the definition of relation, see Unit 19.) In the language of 
section 17.2.1, Unit 17, Logie 11, an equation is an open sentence, which 
becomes a true proposition when an element from the solution set is 
substituted for the variable. 


Thus a differential equation is a relation defined on a set of (real) functions, 
and so its solution set is a subset consisting of functions. For example, one 
element of the solution set of the differential equation 


f+tB=O 
is the function 
fix —Lexp(-2x) (xe R). 
Note that a relation between functions can be written in terms of the 
images under the functions. For example, 
f'+9=0 
is equivalent to 
I(x) + (x) = 0, 
for all x in the domain of f’. We shall refer to the “image form” of a 


differential equation as a differential equation too; in this case it is 
frequently convenient to use the Leibniz notation for the derivative. 


Differential equations arise directly from many basic physical laws and 
are therefore fundamental to the study of considerable parts of science and 
engineering, For example, one of Newton's laws states that 


force applied = mass x acceleration. 


Acceleration is the rate of change of velocity, and velocity depends on 
time. If, for instance, the force applied itself depends on the velocity, as it 
does in the car engine, then we can write the above equation in terms of a 
“position” function and its first and second derived functions, and hence 
obtain a differential equation, Differential equations do not arise solely 
from this kind of background: they also arise from such diverse sources 
as architecture, biology and economics. In this unit we are going to solve 
only a simple form of differential equation involving Q,Q’ and known 
functions. 


Our plan is to concentrate mainly on a few closely related differential 
equations, In section 24.1.1 we shall do the modelling; in other words, we 
shall develop some differential equations from plausible physical situations. 
In section 24.1.2 we shall discuss some basic ideas concerning the solution 
of differential equations. We:shall then investigate graphical solutions of 
differential equations, since the pictorial approach will give us a good 
idea of how solutions behave. This approach will also enable us to appreci- 
ate the numerical approach to differential equations which is introduced in 
section 24.2.4. In sections 24.2.2 and 3 we introduce some which might be 
called formula methods of solution of differential equations. At this 
particular stage the objective of a formula method is to rearrange the 
differential equation in such a way that we can use our powerful tool — 
integration. 

This is the approach to problem-solving which Polya stresses. We ask: “Is 
there a problem like it which we have already solved?” The answer may be 
“Yes” if we do a little rearranging first. 
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24.1.1 The Population Growth Problem 


In this section we take up a problem of population growth considered 
previously in section 7.4 of Unit 7, Sequences and Limits I, and look at the 
increasingly sophisticated mathematical models we can design. 


Let g = Q(t) represent (at any time t € R) the population (in thousands) ofa 
particular species, which we shall take to be human beings, but could 
equally well be viruses, locusts, fish or birds. 


Since q is measured in thousands, it can really take only certain rational 
values corresponding to a whole number of human beings, e.g. 285.632. 
This brings us to our first modification in the process of formulating our 
mathematical model. To talk about growth, or rate of change, we want to 
talk about the derived function, and the derived function was certainly 
not defined for functions with this subset of the rationals as codomain 
since the limiting procedure we adopted in the definition then becomes 
meaningless. This looks like a full stop: but what we do is to take the co- 
domain of Q to be R* in order to make our mathematical model. That is, 
we assume that the variable g can take any positive real values to enable 
us to use the powerful tools of calculus. 


To determine how the population will change with time we need to intro- 
duce functions to represent the number of births and deaths per year. We 
shall call these B and M respectively. The number of births (deaths) 
usually depends on the size of the population, Q(:), which itself depends 
on the time at which it is measured. Thus the composite function B = Q 
will tell us how the number of births per year depends on the time. 


The domain of each of the functions B and M is the subset of the reals 
which represents the population in thousands, and the codomain is a sub- 
set of the reals which represents the number of births (deaths) per year. We 
shall measure ¢ in years. The rate of growth of population is then given by 


DQ=BeQ-MeQ 
or 


Q(t) = B(Q(t)) — M(Q(t)) = Big) — Mig) 
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wd 
Discussion 


The simplest assumption to make is that the number of births and the 
number of deaths per year are both constant, ie. 


B:q*— by (qeR>), 

M:q*— mo (qeR*), 
where by and mo are known numbers. 
Thus we get 

Q(t) = by — mo = ko, say 
It is not our purpose to solve differential equations in this section, although 
this one has obvious solutions. All we are illustrating here is how differen- 
tial equations arise, It is unrealistic to imagine that the number of births 
per year will remain constant; for example, if the population doubles it is 
likely that the number of births also doubles. To translate this into mathe- 
matical terms is to require that the function B be of the form 

B:g—> bq @ER*), 
where b, is a positive number. That is, the number of births is proportional 
to the number in the population. (The number b, is the birth rate in the 
usual official sense; that is, the number of births per thousand per year.) 
Similarly, it is reasonable to assume that the number of deaths is propor- 
tional to the population, i, 


M:q*>—> mq (qeR*), 

where m, is a positive number, The differential equation is now 
Qt) = b, Qt) — m, QW) = k, QO) 

or 


dq 
dt 


where ky = b, — my. 


=kq, 


Now we introduce a further refinement to our model. As the population 
increases, so the available food supplies may become depleted. With a 
greater possibility of disease and disasters claiming a higher number of 
deaths, we would expect that the death rate would increase with increasing 
population 
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We could perhaps represent this mathematically by 
M:gr—>mq + mq? (ge R*) 


where m, is a positive number. For instance, if we take m, = 20 and m, = 
4, we would have 


Miqi— 20g +4q? (qe R"); 
the graph of M is depicted in the figure. 


+ Mia) 


This would imply that with a population of 1 000, i.e. q = 1, the death rate 
would be 24 per thousand per year, whereas if the same population were 
to increase to 10000, ic. g = 10, the death rate would increase to 
200 + 400 
10 
now be 


Ot) = kOe) — m(Q()? 


= 60 per thousand per year. The differential equation would 


or 


dg _, 2 
ae 7 ha mag’, 


where k, = by — my. 


From now on in this text we shall write the “image form” of a differential 
equation in terms of either Q(t) or g, rather than continue to write both 
versions. The first form, explicitly in terms of the images, is useful, since it 
tells us that we are looking for a function as a solution rather than a num- 
ber. The second form, using the Leibniz form of the derivative, is convenient 
because it is more concise. We shall use either form depending on which we 
feel is more appropriate in the context. 


It is interesting to notice that we can obtain some information from a 
differential equation even without solving it. We can write our differential 
equation in the form 


dq _ ky 
ae? mala ‘) 
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ag 
dt 
would decrease until the model became inapplicable or the population 


became zo) we see that, provided 


Assuming k, > 0 {if it were not, — would be negative and the population 


“ is positive, implying increasing population. (Remember that q = Q(t), 


k 

ite. q is a number.) If q = = then the population is static, and if q > a, 
2 

the derivative is negative and the population is decreasing. This tells us 


that Ms is the only stable population. 
2 


This leads us to the final refinement we propose to make here. This stable 
population which can be supported is likely to increase with time, due to 
improvements in medicine, agricultural technology, etc. Let us suppose 
that this stable population increases linearly with time and that we replace 


fe by (kat + k,), where k, and k, are positive numbers. The differential 
2 


equation becomes 


4. m3q(k3t + ks — q) 

Notice that the right-hand side now involves the two variables q and f, 
These variables are customarily given special names; t, the domain vari- 
able of the function Q we wish to determine is called the independent 
variable ; is the dependent variable, since itis the variable in the codomain 
of the function Q we wish to determine. 


Notice also that we have made various suppositions without any real 
justification. What we have is a very tentative model which, on analysis, 
may or may not fit the facts. In a real situation we would need to test its 
consequences against known data before we could attempt to make any 
predictions from it, 


The various differential equations in this introduction are all differential 
equations of the first order. Order is determined by the highest derivative 


< . d" 
present. In other words, if a differential equation contained a or Q"(1) 


and no higher derivatives, we would say it was nth order or of order n, For 
example, the equation 


dq) 2 
“3 ae 


is third order. 
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24.1.2 Basic Ideas about Solutions of Differential Equations 
Consider the familiar problem of solving the quadratic equation 
P-—St+4=0  (reR). 
The solution set of this equation, namely 
{t:? — 5t+4=0,teR}, 
contains two members and is 
{4, 1}. 


Other solution sets in R for quadratic equations may contain two, one (in 
the case of two equal roots) or no members. 


Exercise | 
Write down examples of three quadratic equations whose solution sets 
contain two, one and no real members respectively. a 
Exercise 2 


Describe, in an explicit form, the members of the solution set of the equa- 
tion 


sint=0 (teR) 


Can you count how many members there are? a 


The solution set of a differential equation may be written in a similar way 
to that of an ordinary equation. Consider the differential equation 


Q(=ko (te R) 
or 
Q(t) — ky =0 


The set of all functions which satisfy this differential equation can be 
denoted by 


{0:0 —ky=0 (re R)}. 


In this case each solution to the differential equation is simply a primitive 
function (indefinite integral) of t—> ko, that is 


tr— kot +e — (te R); 


¢ is a constant of integration (see Unit /3, Integration 11). Thus another 
(explicit) form of the solution set is 


{0:0 = t-—+ kot +. ¢ (re R), cE R}, 


which is often more concisely written in terms of images as 


{Q:O(t) = kot + ¢ (te R), ce R}. 
This set can also be written as a relation between variables: 
{Q:q=kot+e (te R),ce R}. 


We see again here the important difference between the solutions of ordin- 
ary equations defined on sets of numbers and the solutions of differential 
equations which are defined on sets of functions. In the former case the 
solution set has real numbers as members whilst in the latter case it has 
(real) functions as members. We can observe this visually by considering 
the graphical representation of the solution set of the equation 


sint =0 (teR), 
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M12 
Discussion 


Exercise 1 
(2 minutes) 


Exercise 2 
(2 minutes) 


(continued on page 8) 
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Solution 1 Solution 1 


For example, 
(i) 7 — 3¢ + 2 = 0 whose solution set is {2,1}. 


(i = 2p whose solution set is {1}. 
(iii) * — 2 + 2 = 0 whose solution set has no real members. a 
Solution 2 Solution 2 


The solution set can be written in the form 
{nx:ne Z}, 


since the sine of any integer multiple of x is zero. There are as many solu- 
tions as there are positive integers. The solutions can be put in one-one 
correspondence with the set of all positive integers as follows: 


1 
n—>2 
—1x-— 3 
2n+—4 
—2n—5 
and generally, 
kn —+ 2k for k > 0, 


kn+—> -2k +1 fork <0. a 


(continued from page 7) 
and comparing it with the graphical representation of the solution set of 
the differential equation 
Q()-1=0 (eR) 
We know from Exercise 2 that the solution set of 
sint=0 (teR) 
is the set of numbers 
{nine Z). 


Graphically, this can be represented as the intersection of the straight line 
which is the pictorial graph of 


t—0 
and the curve which is the pictorial graph of 
t-—+sin t. 


The points of intersection have co-ordinates (nz, 0), ne Z. 


The solution set of the differential equation 
Q(t) -—1=0 (teR) 


may be represented graphically by the family of curves, or set of solution 
curves, given by the pictorial graphs of the functions 


Qit—r+e (teR)ceR. 


The fact that there is a family arises from all the different values that the 
constant of integration c may take. We illustrate the cases c = 2,1, 
0:06, —1,:—2, 


thet-1 


tret-2 


The whole plane is filled by curves (straight lines in this case). Through any 
point, one, and only one, curve of the family passes. Picking out one point 
in the plane picks out one particular curve, which is equivalent to picking 
out a particular constant of integration. For example, consider a popula- 
tion of fish which starts (t = 0) with a population of 2 000 and increases 
at the constant rate of 1 000 per year (difference between number of births 
and number of deaths). 


teet-1 


Then, using units of 1 000, we have an initial point (0, 2) (marked with an X 
on the diagram), and the appropriate solution curve (marked in red on the 
diagram) is the graph of 


Q:t-—+r+2 = (teR) 
since our general solution (i.e. Q:t—> t + c) implies 
O00) =0+e 


2=0+e. 
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In a differential equation obtained as a mathematical model of a real 
situation, we usually know an initial assigned value, as in this fishy 
example. Here, if the modelling is satisfactory, the one particular solution 
curve through the point 


(initial time, initial population), 


can be used to determine the future population. This process of picking a 
particular curve by the known initial condition is an important feature of 
the mathematical modelling in this situation. 


Exercise 3 Exercise 3 
- (3 minutes) 
Sketch a typical subset of the family of curves which illustrates the solution 


set of the differential equation 
Q(t) —-2r=0 (te R). 


Indicate the particular solution curve which satisfies the initial condition 
Q(1) = 2. What is the numerical value of the constant of integration then? 


a 

Exercise 4 Exercise 4 

: (2 minutes) 
Rewrite 

Q(t) = 2 

in a form involving 
(i) functions only, 
(ii) variables only. a 


A differential equation, like an ordinary equation, may have no members Discussion 
in its solution set. Thus 


(eUor + (for +1=0 — (reR)} 
is an empty set, just as 
{i -26+2=0 (te R)} 


is empty. We have just seen that the number of solutions of Q(t) — 2 = 0 
is not finite. Also a differential equation can have a finite number of func- 
tions in its solution set. Thus 


{£10 + (Fo? = 0 (te R)} 
has one member, 

Yi en hd (re R), 
just as 

{tz —2r+1=0 (te R)} 


has only one member. 


Summary ‘Summary 


Solutions of differential equations are functions. The pictorial graphs of 
these solutions form a family of solution curves. In our example we used 
one initial condition (one point in the plane), to pick out one particular 
function (one particular solution curve). 


24.2 SOME METHODS OF SOLVING DIFFERENTIAL 

EQUATIONS 
24.2.0 Introduction 
In this section we look at three ways of finding solutions to first order 
differential equations. Each has its advantages and disadvantages and 
we shall point out some of these as we go along. 
In section 24.2.1 we discuss a graphical method of solution and give an 
example. 
In sections 24.2.2, 3 we pick out from the repertoire of formula methods of 
solution two which will give exact solutions to first order differential 
equations of certain types. (The solutions are exact in the sense that we 
can specify precisely the set of functions which form the solution set, 
rather than give them in tabulated form as with a numerical method.) The 
usefulness of formula methods lies in the fact that they produce a more 
“general” solution than the graphical method, just as the solutions of the 
quadratic equation 

ax? + bx +0¢=0 
are written more generally as 

a —b + \/b* — 4ac 

< 2a ' 

than by saying that if a = 1,b = —2 and ¢ = 1, then x = 1, because the 
general solution can be evaluated for many different numerical values of 
a, b,c. 
Why discuss only two formula methods? The reason is that the two we have 
chosen are simple enough to grasp and learn to use and are representative 
of the exact methods available. It is not our intention to list all the possible 
“recipes” for solving equations —after all, we do not intend that this 
should be a cookery book. Rather, we hope that you will gain a general 
idea of what is involved in finding exact solutions of differential equations, 
and will be able and prepared to learn other methods if and when you need 
them. 
All the exact methods necessarily rely on the methods of integrating and 
differentiating functions we discussed in Units /2 and /3. We shall quote 
them explicitly as required. 
Finally, in section 24,2.4 we give one example of a numerical method of 
solution. 


24.2.1 Graphical Methods 
We wrote down the solution of the equation 

sint=0 
straight away, since we knew that 

sin(nn)=0 (ne Z), 
and no other values of r in the domain R of the sine function are mapped 
to zero. But suppose that we have a less tractable equation of the form 
f(t) = 0 to solve; for example, the equation we had for the “omelette 
problem” on page 34 of Unit 2, Errors and Accuracy: 


‘ 2a 
Pint = 0. 
There we used graphical methods. We calculated the images of the func- 
tion 


E 2n 
ed Se (teR) 
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(continued on page 13) 


Solution 24.1.23 


4a(t) 


< 
is 


The solution set is 
{2:0 =r +e (te R)ceR}. 
Using the initial condition 
Ol)=2=1+e 
we obtain 
c= 
and the particular solution curve corresponds to 
An=r+. 


The family of curves is, in this case, a family of parabolas. All the curves 
look similar in the same way that a biological family may have similar 
traits, That is one reason why it is convenient to use the word family to 
describe this particular collection of curves. a 


Solution 24.1.2.4 


(i) Q=t-—o2 (teR) 
or DQ=t-—>2t (teR) 
a dq_, 
(ii) a= a 
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‘Solution 24.1.2.3 


Solution 24.1.2.4 


at a few selected points in the domain, joined them up smoothly (on the 
correct assumption that the function involved is continuous) and the 
value of t where the sketched curve crossed the x-axis was the approximate 
solution we wanted. 


tt-sint- 2 


(3.14,1.05) 


“+ 


(2.09, -0.87) 


(1.05, -1.91) 
),-2.09) 


Rough sketch of t++t~sint-2Z 


To improve the accuracy, we effectively magnified the portion of the graph 
in the neighbourhood of the approximate solution, by calculating new 
images of the function and drawing that part of the graph more precisely. 


We can adopt a similar approach for the solution of first order differential 
equations; this approach will be useful for those that we cannot solve by 
the standard methods that we discuss in this text. The last equation in the 
population growth example of section 24.1.1: 


di 
7 = mygk3t + ks — q) 
is one of this type. We shall however use the previous equation: 
dq 
ae 7 B14 — maa? 


to illustrate the graphical method, Of course, to solve it graphically we 
must use actual numbers for k, and m3, This illustrates the restrictive 
aspect of the graphical or numerical approach — not only do we have 
to specialize the problem in this way, but, in any case, we do not get a 
general expression for the solution. Suppose we choose k, = 2andm, = 1. 
The differential equation is then 

dq 

dt 
and we want to find the function, or functions, Q, such that q = Q(i), 
which satisfy this equation. We shall be content to have such a function, or 


functions, specified graphically. We shall call the corresponding curves 
solution curves. 


=2g-¢ 


d 
From Unit /2, Differentiation I we know that 4 = Q(t) represents the 
slope of the graph of the function Q, that is, the direction in which the 
solution curve is pointing at any point (r, q). For our differential equation 
we can calculate the slope ofa solution curve at any point ; it is always given 
by 
2q —q?. 

Thus at the point (0, 1), representing a population of 1000 at time 1 = 0, 
the slope of the particular solution curve is 


2x1-P=1 
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(continued from page 11) 


Since 2q — q? is not explicitly dependent on ¢, the slope at any point (t, 1) 
is always 1. In terms of the population growth which the differential 
equation represents, this means, as we would expect, that the rate of 
increase is explicitly dependent on the size of the population and not the 
time at which we are measuring it. In terms of the solution curves, this 
means that for q > 0 we have a solution curve through every point (t, 4), 
and foreach q the tangents to the solution curves at(t, g)are parallel forall r, 


Let us plot slopes at a selection of points on an appropriate graph. Remem- 
ber that the slope is the ratio: 


vertical distance 
horizontal distance’ 


Let us therefore fix a horizontal distance of, say, 4 small squares on the 
graph paper, and then draw the arrow for a particular point to indicate the 
direction of the solution curve there. For example, for the point (1, 0.5) we 
need to draw an arrow to represent a slope of 0.75. The diagram indicates 
how this is done. 


_(1,0.5)-——4 pel 


+ 


The actual numerical values of the slopes at some more selected points are 
displayed in the table below. As we have noted, the value of the slope does 
not depend on the value of t; so all the columns are the same. We have 
nevertheless set out a table, both to give you an idea of what happens in a 
more general case and to produce a record which is similar to the graphical 
picture, 


ane | 
3 -3 -3 | -3 -3 -3 -3 
2.5 - 12s | =1.25 —1,25 P= 1.25 —1.25 —1.25 
z 0 0 0 0 0 0 
| | 
1.5 0.75 0.75 0.75 0.75 0.75 0.75 
1 1 1 1 1 1 | 
=f 
0.5 0.75 0.75 0.75 0.75 0.75 0.75 
0 0 0 0 iF 0 i 0 [ 0 
q 
0 0.5 1 LS 2 Pe 
t 
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If we draw the corresponding arrows, we obtain the following diagram: 


There are several points to note about the table and the diagram: 


(i) Since all the numbers in any one row are the same, wherever solution 
curves cross a given horizontal line they must all be pointing in the 
same direction, 

(ii) Normally we are interested only in that portion of the plane for which 
q > Oand ¢ > 0 (called the first quadrant) since this is the portion 
which is meaningful in the context of this problem. 

(iii) We could go on to calculate the slope at more and more points in 
the first quadrant and make the set of arrows as numerous and dense 
as we please (physical width and length of pencil marks allowing). 

(iv) In the same sense that we sketch a curve through a few points, so too 
we can sketch a curve, or curves, using the arrows as guides. The 
following diagram shows a few typical solution curves. All we have 
to do to sketch a solution curve is to start from a given point and 
make sure that the direction of the curve is roughly parallel to all the 
arrows it passes close to. 


(v) If we wish to improve the accuracy, we must use more arrows, just 
as we would use more points to draw an ordinary graph. 

(vi) What we have produced is an approximation to a selection of curves 
from the family of curves which represent the solution to the differen- 
tial equation. We have thus, albeit laboriously, got the picture of what 
the family of solution curves look like. If you look at the set of arrows 
at a distance you get the impression of how the pattern “flows”. 

(vii) For our particular differential equation, we can see by substitution 
that g = 0 is a trivial solution: if the population is zero initially, 
mathematical and other reasoning tells us that the population will 
never be anything but zero. When q > 0, all the solution curves 
appear to approach the same numerical value 2 as 1 becomes large. 
20000 is thus the stable population figure. If there were more than 
2.000 people originally they will reduce to this number: if less they 
will increase to this number. As we mentioned in the last section, 
we could have deduced this particular piece of information without 
solving the differential equation. For the derivative, Q’(r), is zero when 

2q-4q? =0, 
thatis when q = O(and we have dispensed with that)and when q = 2. 
The value q = 2 is never attained in our present model except in the 
special case when the population starts with 2 000. 

(viii) To select the particular solution curve representing a particular case 
of interest, we need more information; e.g. an initial condition. 
Suppose we start with 750 in the population; that is, q = 0.75, and 
choose to measure time from that instant, so that g = 0.75 when 
t = 0, We start to draw our curve from this point and produce the 
approximate solution curve indicated by a red dotted line in the 
previous diagram. After 2 years* we see (from our curve) that the 
population of 2000 has almost been attained. After one year we 
estimate that there would be | 700 in the population. 


The graphical method works for all first order differential equations which 
we can put in the form 


Agi. 
aes S(t, q) 


with certain conditions on /. All we have to do is to evaluate the images 
under f at appropriate points in the plane. The conditions on f which 
ensure that a solution exists are too complicated for the Foundation 
Course and will be dealt with in a later course. 

The graphical method can be a very useful method in cases which cannot 
be solved by any other means. It can be inaccurate (as can curve-sketching 
of any type) and it is certainly laborious, but one advantage is that we get a 
good qualitative idea of the shapes of all the solution curves, that is, the 
solution set as a whole. With this “panoramic view” we can choose 
particular points of interest and find numerically the particular solution 
curves which pass through them. Another advantage of the graphical 
method is that it gives an introduction to numerical methods of solution of 
differential equations. We look at these in section 24.2.4. 


Summary 
In the graphical approach we use the geometric information given by the 
differential equation 

dq 

de Lh 


namely the slopes of the solution curves at all the points (t, q) of the domain 
of f. to sketch the family of solution curves. 


* Better late than never! 
Were you quicker than the author? 
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24.2.2 Formula Method 1: Separation of Variables 
This method is useful only when we can “disentangle” the variables in the 
equation. 
We illustrate the method by an example. Consider the equation 

QO) = Q(t) (rE R). 
We regard t and Q(r) = q as variables, and we can “disentangle” them by 
rearranging the equation in the form 

ow _, 

Or) 
provided that we exclude any ¢ such that Q(t) = 0 from the domain of Q. 
Writing this in function (rather than image) form we get 


1 
~x DQ=Ir— 1 (te R). 
@ {4} ) 


4 1 2 3 
Since = x DQ and t-—+1 are equal functions, we can equate their 


Q 
primitives provided that we choose the constants of integration correctly. 
We therefore try to “integrate both sides”. 


We can integrate the right-hand side to get the set of primitive functions of 
e 
2 
be so obvious.* In Unit /3, Integration I] (section 13.2.4) we obtained the 
following formula for the primitive function of a composite function: 


fiek = Dk ({s}>« 


Comparing this with our equation, we see that we need to choose k = Q 


the form 1+—+ — + e, ce R, but the integral of the left-hand side may not 


1 
and go k =~, that is 
ee 


u— 
‘J t 


Then from our table of integrals we have 
fe=roine+e, (teR*) 
and 


({s)- =t-— In Olt) + cy; (Qe R*), 


where c, is any real number. Notice how another condition has crept in 
here, that is, Q(t) > 0 (because the logarithm function has domain R*). We 
could deal with Q(t) < Oifit were ofany interest, but we would have to do it 
separately. When using standard results, as we are here, we must always 
be careful to use them in the appropriate circumstances. 


Putting the bits together, we have 
2 


t 
t— In Q(t) + ¢, = Se hie 


This specifies the set of functions Q which belong to the solution set of the 
original equation, but in an inconvenient form: the dependent variable 
q = Q(t) is not expressed explicitly in terms of the independent variable t. 


“In fact, we are on familiar ground ; see Exercise 2, section 12.3.2 of Unit 12. 
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2422 
en est 


So we try to reorganize. First we use the simpler form in terms of images: 


2 
Ing+¢, a Gee 


e 
Ing=> + (¢ — ¢) 


Since c and c, are any real numbers, ¢ — c, is any real number. We can, 
therefore, either write de R for ¢ — c, or just drop the c, . (In general, we 
do not need to introduce the constant of integration on both sides of an 
equation, since we know that if two functions are equal, their primitives 
differ by a constant.) So now we have 

2 


t 
Inqg=z+d 


Remembering that In is a one-one function and exp is its inverse, we have 
P 
=exp(~t+d 
nike °(5 


We now have the dependent variable q = Q(t) expressed explicitly in terms 
of the independent variable f, and we could say that we have a satisfactory 
solution to our problem. In fact, we can tidy things up a little further if we 
remember that 


exp x = e* 


So 


ao (5 + 4) = te x A 


If we write A = e*, the general solution becomes 
q = Ae”? or A exp (5) 


where A € R”, since exp d cannot be negative or zero. 


The family of solution curves is illustrated below. 


4ait) 


A=2 
A=1 
A=0.3 


> 
t 


Graphs of Q(t)}=Aexe( ) 
for A=2,1and0.3 
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We have solved this equation in small steps, and it is worth noticing that 
we have sprinkled conditions on the way. Originally we required Q(t) # 0, 
and then we required Q(t) > 0, which superseded the first condition. 
Finally, we have 


Fe) 
Olt) = Aexp (5] 


and since Ae R*, the condition Q(t) > 0 is automatically satisfied. Also 
notice how the constants of integration, cand ¢, , which originally appeared 
added on to the ends of our expressions and were unrestricted, became 
“entangled” as a result of our manipulations, so that our final constant, A, 
is qualified (A > 0). With practice, the formal manipulation can become 
deceptively easy, but the most adept manipulators often forget to consider 
the consequences of what they are doing. 


And yet, having been so ponderously careful, what if we had not been so 
careful? You can easily check that 


? 
" 
with AéR~, also belongs to the solution set of the original differential 
equation. (The solution curves are similar to those described in the last 
figure but they are all upsidedown below the t-axis.) So perhaps a little 
carelessness pays. Well, it may do, just as in any other scientific field ; but it 
must be disciplined carelessness. In this case, the discipline lies in checking 


that our formal manipulations, without conditions, have not caused us to 
include in our solution set functions which are not in fact solutions. 


Olt) = Aexp 


To generalize the approach using the separation of variables, we need to 
be able to recast the original differential equation in the form 


(g°Q) x DQ =h, 
or, in terms of variables (images); 
dq _ 
807, = hi), 


and then we hope to be able to integrate. Herein lies the reason for the 
name of the method ; we “separate” the q’s to one side of the equation to 
enable us to integrate. This means that the original differential equation 


must be expressible in the form 


ag RO: 
dt gq) 


Exercise 1 
There are some fairly sophisticated conditions on g and h for this method 
to work, which we shall not discuss in this text; but what obvious re- 
striction must we place on g? a 
Returning to the form in Equation (1), we have 

(g2Q) x DO=h 
Using the result for integrating a composite function quoted above, we 
get 


fe-@x ve= [r+o—o 


‘Equation (1) 


Exercise 1 
(2 minutes) 


Main Text 


(continued on page 20) 


Solution 1 


g(q) may not be zero, ic. we must exclude those values of ¢ for which 


(Q(t) = 0. a 


(continued from page 19) 


that is 


({s]-o-[a+u—o 


Thus, if we can obtain f gand f h from appropriate tables (perhaps after 


further manipulation), we can solve the differential equation. 


Summary 
If 
h 
DQ = —. 
2 g°Q 


then the solutions satisfy the equation 


[/« 


The Leibniz form of the separation of variables rule is outlined below. 
If 


(g(Q(t)) 4 0. R). 


2a=[h+e—o ceR. 


dq hit) 
> =— #0), 
i ae (g(q) # 9), 


then the solutions satisfy the equation 


J etoraa = fino at +e 


This is a form you will probably find in textbooks ; you may find this form 
easier to remember. 


Exercise 2 


Find the function Q which is a solution of 


di 3 
Cr =2t (teR) 
where q = Q(t), and which satisfies the initial condition 
Q(0) = 3 1] 
Exercise 3 


Find the solution set of 


Eel 
dt qt+a 


where g = Q(t), and a and k are positive numbers. State any conditions 
you find necessary in the manipulations. a 


Exercise 2 
(3 minutes) 


Exercise 3 
(3 minutes) 


FM 24.2.3 


24.2.3. Formula Method 2: Integrating Factor 24.23 


We can now find formula solutions for first order differential equations Main Text 
of two types ‘a 


(i) OH) = gir). 
by straightforward integration, and 


a A 
(ii) Q() = now h(Q(t)) # 0, 


by separating the variables, provided the resulting integrals can be found 
with reasonable ease. Obviously there will be cases where the equation 
does not admit of the particular form required in the separation of 
variables method. For example, 


Q() = —Qlt) + Pexp(—1) 


cannot be put in the required form. We can separate the variables in the 
sense that we can get all the dependent variables on one side, i.e. 


Qt) + Olt) = P exp(—n, 


but we cannot get the left-hand side into the form (g° Q) x DQ. However, 
we can adopt the same idea. Our separation of variables method relied on 
our being able to use known results about the integration of composite 
functions, and these were themselves obtained from our result for differ- 
entiating composite functions. Now we also know how to differentiate 
the sum and product of functions. The “sum” result is so obvious that it is 
unlikely to give us any technique for solving differential equations which is 
not itself obvious. But the “product” result (Unit /2, Differentiation 1, 
section 12.2.4) 


Dif * g) = (Df) x g + f x (Dg) 


is not so obvious. Suppose that we assume that the left-hand side of our 
differential equation is of the form D(f x g), for some choice of fand g. 
Then 


(Df) x g + f x (Dg)=Q'+Q 
Now we can choose f = Q; then we have 
Qxg+Qx (Dg)=O'+O 
i.e. in image form 
QMie(git) — 1) + Q(Dgit) — 1) = 0 
If Q and Q’ are linearly independent,* then 
2x0 + PQ’ =t-—+022=f =0 


2Q(t) + Q(t) = 0x =Pp=0 
In this case, we have no further choice, in the sense that 
Dg =t-— 1 (teR) 
and 
g=t—1 (teR) 


and this is plainly impossible. However, now that we have started on this 
way of thinking, let’s continue with it. We obtained the form 


Q(t) + Olt) = P exp(—) 


by separating out the terms involving the dependent variable q = Q(t); 


* See Unit 22, section 22.2.2. (continued on page 23) 


Solution 24.2.2.2 


The equation is already in separated form. Comparing it with the 
equation 


dg _ He) 
dt g(a) 
we see that 


(a) = 4. h(e) = 2 
Therefore suitable primitive functions of g and h respectively are 


2 
fa *q=q sc and fe — 2 = tr ? 


So the solution set is 
¢ 
foarte werveer}. 


The particular value of ¢ determined by the initial condition is given by 
4(Q(0))? = 07 + ¢, where Q(0) = 3, 


ie. 
c=¥ 
Therefore the particular solution we are interested in is given by 
(OP _ 2,9 
has ek 


But this determines two possible functions Q: 

(i) Q:t-— 2° +9 

(ii) Q:t-— — /2P +9 

Remembering again that Q(0) = 3, Le. Q(0) > 0, we see that the required 
solution is 


Q:t-—> /20? +9 a 


Solution 24.2.2.3 


A first condition is obviously that q + a # 0, i.e. Q(t) # —a. In separated 
form the equation becomes 


qtadq 


es =k: 


q dt 
and we have to impose the further condition q # 0, i.e. Q(t) # 0. Referring 
back to the summary, we have 
qt+a a 
) =—— = 1 +,-, i) = —k 
aq a = 


Suitable primitive functions are 


fe-a—a+aing (assuming q€ R*), 


h=t-——kt 
Therefore, the solution set is 
{0:0() +alnQ(t)=—kr+e = (te R), cER, OER } 


Notice that Q(t)¢R* takes care of the previous restrictions on Olt). 
This time we cannot get the solution set in an explicit form, that is, we 
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‘Solution 24,2.2.3 


cannot get q, the dependent variable, explicitly in terms of t, the inde- 
pendent variable, and thus it is difficult to determine Q(t) for a given t. 
In fact, we could use one of the iterative methods we introduced in Unit 2, 
Errors and Accuracy, or evaluate Q(t) at a selection of values of t and then 
use interpolation to find other values of Q(t). On such occasions it may be 
just as simple to solve the differential equation by a numerical method 
straight away. a 


(continued from page 21) 


and then we tried to see if the left-hand side had product form. The latter 
part is the important bit, so let’s concentrate on that. Can we find a left- 
hand side, variables separated or not, such that we can recognize a 
product form? Assuming that one function in the product is Q, say f as 
above, then we require a form 

Ox g+Qx Dg. 
Comparing this with what we have actually got in the differential equation, 
we must somehow find a way of multiplying Q’ by a function g, and Q by 
its derivative Dg. Now we can’t multiply Q’ by one function and Q by 
another without changing the differential equation: we must multiply 
them both (and of course the right-hand side) by the same function. This 
function must therefore satisfy the equation 

Dg=8 
Now any function g satisfying this equation will do, and we know one 
solution, the exponential function, since 

D(t-—+ exp t) = t-—rexpt 
So we multiply both sides of our equation by exp ¢ (which luckily is not 
zero anywhere, so we are not likely to produce nonsense). We obtain 


(exp NO) + (exp NQ(t) = (exp d(exp(—0) = 


Diexp x Q) =1-— P, 


or 


rid 
SPQ OS Fs, 


from which we can get Q or Q(t) = q explicitly as 
3 
q = exp(—t) x (5+ (| 


So the solution set is 
3 


2:01) = exp(— x [+e] (eRycert. 
3 


Before we discuss some general aspects of this method we suggest you try 
it for yourself in the following exercise. 
Exercise | 


Find the solution set of the equation 


C= a 
Pan te (teR,t # —a) 


where a is a positive number and q = Q(t). a 
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Exercise 1 
(4 minutes) 


Solution | 


As before, we assume that we can get the left-hand side of the equation 
in the form of a derivative of a product, i.e. we assume that it can take the 
form 

digi) 


4a x a0) 4) + 
at * BO) = G8) + aa 


To get the left-hand side of our differential equation into this form, we must 
dg q d(git)) 

ty —1_ py 7 

multiply r by a(t) and oar by ze x (1 + a). These two must be the 


same, i.e. we need a function g such that 


digit) _ 
+ a)—7 = 

We can separate the variables in this equation in the sense of our first 
method, to obtain 

1 digit) _ 

g(t) dt t+a 
Weare looking for any one solution of this equation so we take the constant 
of integration to be zero. 
(We have seen the left-hand side before.) 
Integrating, we get one solution 


In(g() =In(@e+ a) (@ +a > 0) 
or 
gi)=t+a 
So multiplying both sides of the original differential equation by t + a 
(which again is not zero, since t + a > 0), we get 
dq 


( + oF, 


+q=t +a) 


(Of course, this rearrangement may have been obvious to you from the 
start: we have given a formal approach just in case it wasn’t.) We have 


d nme: 
qt + a))=t + at, 


so that 


which gives the solution set 


3 2 
Hl Lali cJirert # -aycer}. 
tta 2 


{o :Q(t) = 


(In the process of solution we required ¢ + a > 0, but this condition is not 
essential: it is easy to verify that the above describes the solution set for 
t +a <Oas well.) | | 
The crucial step in both the example in the text and Exercise | is the 
multiplication by a suitable function g, so that the left-hand side of the 
equation can be written as the derivative of the product of g with Q, ic. 


Dg x Q)=80' + '0 


If we are going to adopt this as the general form, then it implies that we are 
dealing with an equation of the form 
g h 
y+=Q=-, g(t #0, 
2 g 2 g 
except that we are not going to be told what g is. 


24 
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So let’s start with a given equation of the general form 
O+PO=R, 


where P and R are known functions. Comparing this with the previous 
equation, we see that g is determined by 


Sup 


g 


Remembering that we are looking for just one solution, we ignore the 
constant of integration, and get 


ineg=[P g(theR® 
or 
g=expe f ra 
or 
g(t) = exp ({ riya) 
Summarizing, we get the following rule: 
To solve an equation of the form Rule 


Q(t) + P()O(t) = Rie), 


where Pand R are known functions, multiply throughout by the integrating 
Jactor 


g(t) = exp (J P(t) a) 


and integrate directly. 

Let us try an example. We shall find the solution set of the equation 
Q(t) + Q(t) = 4% exp(—1?) (te R) 

To compare with the rule, we take 
Pst——o2t 

and a simple integral is 
et 


Therefore an integrating factor is 
g(t) = exp (J P(t) a) = exp(t?) =e", 
and we have 


(eo) =e" Q(t) + 2te" Q(t) (derivative of a product) 


= Q() + 2) 
=e"(4re"*) (from given differential equation) 
= 4 


(Note that the left-hand side is 4 Gntegrating factor x Q(0).) 
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The integral is then 
Ate* =t* +c, 
and the solution set may be written as 
(0:0) =e +0) (te R)cER} 


Exercise 2 Exercise 2 
(4 minutes) 
Indicate which of the following differential equations can be solved 
(assuming that the appropriate primitive functions can be found for the 
integrals involved), by 
A: both the method of separation of variables and the integrating factor 
method; 
B: the method of separation of variables and not the integrating factor 
method ; 
C: the integrating factor method and not the method of separation of 
variables ; 
D: neither the integrating factor method nor the method of separation of 
variables. 


de 

‘) oF = in(ea) (teR*,q = QWeER") 
F dq_ 

(ii) (exp = 4 

(iii) OM = 10) (He R*, ER”) 


(iv) a“ = cos q + cost 


dq_ g+i 
a Pel 
Solve two of the above equations for which you can find the primitive 
functions for the integrals involved. a 


The remainder of this section can be omitted if you are short of time. 

As a last example of a formula method, we solve the population growth Example 
equation (see page 5): 

Oo kg — maa? (teR,qeR*) 

where q = Q(t), and k, and m, are positive numbers. 


Neither the integrating factor method, nor the method of separation of 
variables is directly appropriate, but there are many other methods 
available. 


We shall replace the function Q by the function U, where 


1 
A = Ta 
16, 
ai 
aa 


where u = U(t). 


If you are familiar with the Leibniz notation, you can deduce from the 
last equation that 


26 


FM 24.23 


(Alternatively, using the notation of Unit 12, Differentation I, you could 
I 
deduce the same result from Q = (»—4) eU) 


aq 


Substituting for at 


and q in the differential equation, we get 


eansel eie ae 
wd 


that is, 

7 o —kyu + m, 
or 

du 


at kwam, 


This equation can now be solved by the integrating factor method. An 
integrating factor is 


g(t) = exp (J a) = exp(k,t), 
so the differential equation may be written: 
Sex (ky tu) = my exp (ky) 
and after integrating we obtain 
exp (k,Qu = Fexpkit) +e 
or 
w= + cexp(—kit) (te R) 


Since qe R*, we R*, so that the constant ¢ cannot be negative. Having 
1 , 
solved the equation, we replace u by 5 to obtain the result in terms of our 


original variable q. 


1 om 

—=—— + cexp(—k,t) 

* a5 P 1 

The explicit form of the solution set is thus 
1 
Q:Q(t) = ————__ (te R§),c > 01. 

Ma 4 ce-kt 
ky 


In this case we can determine the population at a specified time t by 
reference to exponential tables rather than by drawing a solution curve. 


Summary Summary 


In this and the preceding section we have described two formula methods 
of solving first order differential equations: the “separation of variables” 
and “integrating factor” methods. If you are faced with the need.to solve a 
first order differential equation which is not amenable to one of these 
methods, you can either refer to other formula methods described in 
standard texts (see Bibliography), or use the graphical method or a 
numerical method. The next section is devoted to introducing numerical 
methods. 
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Solution 2 Solution 2 
(i) D 
(ii) B 


Solution set is 


{oo = weR,cer} 


1 
exp(—t)+c¢ 


If c is negative, there will be the problem that Q(t) is not defined for 
one value of ¢ and so there would be a discontinuity in Q. 


(iii) A 
Solution set is 
{Q: Q(t) = exp (3°? + ) (teR*),ceR} 
or 
{Q: Q(t) = aexp (3°) (teR*),aeR*} 
(iv) D 


(v) B 
To find the solution by the method of separation of variables requires 
us to be able to recognize the primitive function 


1 
i—kea— 
i P+ 
We have not formally covered this in the course: it is, in fact, 


t— are tant (te R) 


which is defined to be the inverse function of 


t— tant ry ee! 
2'2L/° 
(See Exercise 1, section 13.2.6 of Unit /3.) The solution set is then 


{Q:are tan Q(t) = arctant + ¢ (te R),ceR}. 


This could be simplified if there were any further interest in the 
problem. a 
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24.2.4 A Numerical Method of Solution 


There are many powerful numerical methods of solving differential 
equations. As with the other methods described earlier in this text, we 
want to give you some idea of just one such numerical method. It is not 
very accurate, but it will give you an idea of the basis of the numerical 
approach. It is a method which, if suitably refined, could well be used to 
solve a differential equation. It is called Euler's Method. The basic idea 
of the method is rather like hunting for pirate’s treasure. 


Go to the solitary tree on the top of the hill. 
Take 20 paces in an easterly direction. 

Then take 15 paces in a north-easterly direction. 
Etc. 


Euler’s method can be used to solve an equation of the form 
O10 = fla), 


with given initial condition Q(to) = qo, where q = Q(t), f is a known 
function of two (real) variables, and ty and qo are known numbers. 


(fo, 4o) is the pair of co-ordinates of the point (the solitary tree on the 
top of the hill) from which the solution curve starts. 


The differential equation tells us the direction in which to set off; we 
evaluate (fo, qo), and that is Q’(to). 


+9 


’ 
' 
: 
| O 

’ 


‘ 
(ta,42) 


Atssain) 


{to,40) 


ne 


In Euler’s method we can decide how far we should go in that direction 
before recalculating a new direction. Intuitively, the shorter the distance 
we go, the more accurate our estimate of the solution curve is likely to be. 
The first step takes us to some new point (r, ,q,) where Q(t,) = q,. 


Again the differential equation will tell us in which direction to go, since 


Q(t) = M.ay) 


is the slope of the solution curve through the point (t,,q,). We can con- 
tinue for as long as we wish. The resulting path consists of straight line 
segments, which we hope are a good approximation to a part of the true 
solution curve through the initial point. 


Let’s see how this method works by applying it to a particular example. 
In doing this we have a choice to make. We could choose an example of a 
differential equation which cannot be solved any other way, or we could 
choose an example of a differential equation for which the solution is 
already known. We have chosen the latter because the former tends to 


Leonhard Euler 
1707-1783 


involve more complicated arithmetic. Also, using the latter for first 
acquaintance with the method, we can see how well we are doing and get 
the feel of what is going on. We shall estimate the solution curve to 


(0) = 200) 
with the initial condition 
Q(0) = 0.2, 


i.e. the required solution curve passes through the point (0, 0.2). We can 
in fact solve this exactly, either by using the “separation of variables” 
method or the “integrating factor” method (or almost from foreknow- 
ledge!): the required solution curve is the pictorial graph of the function 


OQ, :t'—> 0.2 exp (21) 


We shall now work on the following figure, on which the correct solution 
curve has already been drawn in red, and use a numerical method of 
solution. The “solitary tree” we start from is marked by a cross on the 
graph paper at the point (0, 0.2). We now have to choose the number of 
“paces” or “distance in a fixed direction” we are going to take. Unlike 
the treasure hunt, we do not measure the paces along the path but take a 
fixed step length in the t-direction. (This is the same as in the other approxi- 
mation methods we used in Unit 4, Finite Differences, when using finite 
differences to interpolate, and in Unit 9, Integration I when evaluating 
definite integrals using the trapezoidal or Simpson’s rule.) 


9 


at 


We shall use a step length of 0.2. The direction in which we have to go at 
the point (0, 0.2) is obtained from the differential equation. It is along a 
straight line path whose slope is 


Q(0) = 2 x (02) = 04 


We continue along this path until we have covered our step length on the 
t-axis, i.e, until the ¢ value is 0.2. (We draw it by measuring 4 squares up 
from the horizontal line through (0, 0.2) for every 10 squares along it.) 


Suppose that this path we are producing is the pictorial graph ofa function 
V, then 


vio) =02 
V(0.2) = 0.2 + (slope of straight line segment) x (0.2) 
= 0.2 +04 x 0.2 = 0.28 
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To calculate the new direction in which to go from the point (0.2, 0.28), 
we calculate the slope of the solution curve through (0.2, 0.28). Note that 
we are already “off” the actual solution curve: our change in direction at 
this point is to try to keep “close” to this curve. At (0.2, 0.28), 


Q'(0.2) = 2Q(0.2) 
= 2V(0.2) 
=2 x 0.28 = 0.56 
We therefore continue the graph of V by following a straight line segment 
from (0.2, 0.28) with a slope of 0.56 until we reach a t value of 0.4. Then 
V(0.4) = 0.28 + 0.56 x 0.2 
= 0.392 
and we have determined the next “corner” of the graph of V. 
Continuing the procedure we would find 
V(0.6) = 0.392 + Q(0.4) x 0.2 
= 0.392 + 2V(0.4) x 0.2 
= 0.392 + 2 x 0.392 x 0.2 
= 0.549 


Exercise 1 


By using a step length of 0.1 instead of 0.2, determine the first four line 
segments of the estimated corresponding solution curve, call it V,, and 
sketch it on a diagram. Does the graph of V, appear to be “better” than 
the graph of V, in the sense that it is nearer to the graph of the actual 
solution? = 


Exercise 2 


You may be wondering how to determine the estimated solution curve for 
t <o, that is, te R~ in this case. How would you determine V(—0.2)? 


Now let us generalize and suppose that we wish to find an approximate 
solution to 
Q = f(a) 
for t in the interval [a, b], given the initial condition 
Qa) = go 
For simplicity, we divide the total interval into n equal step lengths 


b-a 
n 


h= 


and label the ¢ values 


to=a 
t=ath 
t,=a+rh 
t,=a+nh=b 
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Exercise 1 
{5 minutes) 


Exercise 2 
(5 minutes) 


(continued on page 33) 


FM 2424 
Solution I Solation 1 


The tabular points for V, are approximately (0, 0.2), (0.1, 0.24), (0.2, 0.288), 
(0.3, 0.346), (0.4, 0.415), 


+9 


It appears that the graph of V, is better than the graph of V as an approxi- 
mation to the actual graph. a 


Solution 2 Solution 2 


Use Q’(0) and march “backwards” until ¢ = —0.2, V(—0.4) and so on 
may be determined in a similar way. 


| 
| 


t=-0.2 | t=0.; 
: if 
eg. 
V(—0.2) = (0) — 0.200) 
= 0.2 — (0.2) x (0.4) 
= 0.12 1] 
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If Vis the estimated solution we have (continued from page 31) 
Vito) = do 
V(t) = qo + hQ'(to) 
= do + hf (to, 40) 
Then 
V(t2) = Vit) + h(t) 
=v, + Af(t,, 04), where v, = Vit,) 


and so on. We obtain successive values by repeating this process and this 
can be represented in the following flow diagram. 


ves 


This determines a set of (n + 1) points on the approximate solution curve. 
Notice that this method corresponds to using 


An 
h 
This can be seen as follows: if 
DQt) = f(t, 4) 
is replaced by 


MV) = Sets 


g as an approximation to Dg. 


we get 


V(t, +h) — Vie) 
ADAMO) _ 96,0), 
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that is 
V(t,+1) = Vit.) + Af. Uh 
or 
O41 =O, thf, 
where 
fp =Sltrstd v= Vit). 
Exercise 3 


The differential equation we used to illustrate Euler’s method was 
Q() = 200) 
with the initial condition Q(0) = 0.2. 
For this equation, show that if the step length is h and we start from t = 0, 
then 
v, = 0.2(1 + 2hy 


Suppose we are determining the solution curve in the interval [0, 1) and 
using n lengths; show then that 


2\" 
+3) 
n, 


What is the exact value of Q(1)? What is jim v,? What do you conclude 
large 


about the sequence of estimated solution curves produced as n increases? 
a 


v, = V(L) = 0.2) 


The last exercise brings us to the problem of the accuracy of Euler’s 
method. Any detailed discussion of this would take us too far: it is 
technically very complicated, although we already have some of the 
necessary concepts, such as error interval, error bounds and the correction 
terms in the Taylor approximations (Unit 14, Sequences and Limits 11). 
If you are interested, see the book by Macon listed in the bibliography. 
But, be warned, it goes well beyond anything we have done in this course, 
and will prove very difficult at this stage. 


Summary 

We have seen in this section an outline of one of the most elementary 
numerical methods of solving first order differential equations. There 
are many more advanced methods which increase the accuracy of the 
estimated solution curve. For example, the Euler method predicts the 
next point from information at the previous point. This new knowledge 
(the co-ordinates of the next point) enables us to go back and correct 
our estimate of the slope to use over the previous step length. This is 
just like an iterative process —a first guess, followed by an improvement. 
Such methods are called predictor —corrector methods and lead to 
important improvements in accuracy. 
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Exercise 3 
(5 minutes) 


24.3, CONCLUSION 


When introducing differential equations we are always faced with the 
problem of the mode of presentation. Should we spend a considerable 
amount of space and time demonstrating how the mathematical modelling 
of a wide variety of problems (such as absorption of drugs into the system, 
aircraft flight paths or water seepage) produces differential equations, 
or should we concentrate on the significance of the solutions of the differ- 
ential equations? 


We have tried to compromise by selecting a single problem, the population 
growth problem. We set up the appropriate differential equations and 
then used this problem as a framework in which to discuss the graphical 
method, two formula methods and a numerical method of solution. The 
graphical method enables us to see how solution curves behave, but it 
only rarely has a direct practical use. The formula methods were funda- 
mental for a long time to finding any solutions at all (and still are to a 
considerable extent at an elementary level). Numerical methods in general, 
that is, methods based essentially on the ideas in section 24.2.4, are very 
powerful methods because of their wide application; however, they 
always have associated with them the problem of determining the accuracy 
involved. 


Now for a look ahead. There will be another unit on differential equations 
in which we shall look at an example of a “second order” differential 
equation; that is, a differential equation involving the second derivative 
of the desired function. There we shall model the specific problem of 
vibration. 
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Solution 24.243 
1+ Af (t,-1+%,-1) 


v,-1 +h x 20,-1 


since, in this case, 

S(tp—15 P= 1) = 20,-1 
Therefore 

v, = (1 + 2h)v,-1, 
and by repeated substitution we get 

v, = (1 + 2hYvy = 0.2(1 + 2hy 
For the interval [0, 1), 

ads 
so that 


2\" 
V(1) =v, = o2(1 + °) 
n 


Now for the actual solution Q, 
Q(1) = 0.2 e* 
and 


lim v, = lim oa(1 + ‘} 
niarge nlarge n, 


which, from Unit 7, Sequences and Limits I, is also 0.2 e*. The two values 
are thus the same and we conclude that the limit of the sequence of 
estimated solution curves coincides with the actual solution curve at 
both its ends if n is sufficiently large. This suggests that we can get as 
close to the solution curve as we please in the interval [0, 1), although we 


have only checked the end-points. 
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Solution 24.2.4.3 


NO TEXT 


NO TEXT 


NO TEXT 


NO TEXT 


Title of Text 


Functions 

Errors and Accuracy 
Operations and Morphisms 
Finite Differences 


Inequalities 

Sequences and Limits I 
Computing | 
Integration I 


Logic I— Boolean Algebra 
Differentiation I 
Integration II 

Sequences and Limits II 
Differentiation II 
Probability and Statistics I 
Logie Il — Proof 
Probability and Statistics I] 
Relations 

Computing II 

Probability and Statistics III 
Linear Algebra I 

Linear Algebra II 
Differential Equations I 


Linear Algebra IIL 
Complex Numbers I 
Linear Algebra IV 
Complex Numbers II 
Groups I 

Differential Equations II 


Groups II 

Number Systems 
Topology 

Mathematical Structures 
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